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^ ! 1. Introduction 

0^ 

O ■ Much work has been done on the 2-dimensional Ising model (IM) during the past over 
60 years. Many analytic results for the partition function and correlations have been 
^ ■ obtained. These have greatly contributed to establish our present understanding of 
■ continuous phase transitions in systems with short range interactions [1, 2, 3, 4, 5, 6]. 
Recent overviews with many references are given e.g. in [7, 8, 9]. Many rather different 
mathematical approaches have been used, so that already 30 years ago Baxter and 
Enting published the "399th" solution for the free energy [10], see also [11]. Spin-spin 
correlation functions can be written as Pfaffians of Toeplitz determinants. Most work 
has focussed on the thermodynamic limit and scaling properties since these give contact 
to field theoretical results and to beautiful Painleve properties [5, 6, 12]. 

Only during the last decade more attention has been drawn to correlations and spin 
matrix elements (form factors [13]) in finite-size Ising systems [14, 15, 16]. Nanophysics 
experimental arrangements often deal with systems where the finite size matters. 
Recent theoretical work on the finite-size IM started from Pfaffians and related Clifford 
approaches. In [17] it has been pointed out that one may write completely factorized 
closed expressions for spin matrix elements of finite-size Ising systems. One goal of the 
present paper is to prove the beautiful compact formula conjectured in Eq.(12) of [17], 
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sec (129). For achieving this, we introduce a method which has not yet been apphed to 
the Ising model: Separation of Variables (SoV) for cyclic quantum spin systems. Our 
approach is the adaption to cyclic models of the method introduced by Sklyanin [19, 20] 
and further developed by Kharchev and Lebedev [21, 22]. We also make extensive use 
of the analysis of quantum cyclic systems given in [23]. 

Little is known about state vectors of the 2-dimensional finite-size IM. Only partial 
information about these state vectors can be obtained from the work of [3]. Recently 
Lisovyy [24] found explicit expressions using the Grassmann algebra method. Here we 
shall present our SoV approach [25, 26, 27] which gives exphcit formulas for finite- 
size state vectors too. However, these come in a basis quite different from the one 
used in [24]. We shall calculate spin matrix elements by directly sandwiching the spin 
operator between state vectors. Factorized expressions result if we manage to perform 
the multiple spin summations over the intermediate states. 

The prototype of a general A^-state cyclic spin model is the Baxter-Bazhanov- 
Stroganov model (BBS) [28, 29, 30], also known as the r^^^-model. The standard IM is 
a very special degenerate case of the BBS model. In order to avoid formulating many 
precautions necessary when dealing with the very special IM, we shall develop our version 
of the SoV machinery considering the general BBS model. We chose to do this also 
because of the great interest in the BBS model due to the fact that its transfer matrix 
commutes with the intcgrable Chiral Potts model (CPM) [31, 32] transfer matrix [29, 33]. 
Obtaining state vectors for the CPM is a great actual challenge [34, 35]. Although the 
eigenvectors for the transfer matrix of the BBS model with periodic boundary condition 
are unknown for N > 2, explicit formulas for the eigenvectors of the BBS model with 
open and fixed boundary conditions have been found [36, 37]. 

This paper is organized as follows: In section 2 we define the BBS model and 
its Ising specializations. In section 3 we discuss the Sklyanin SoV method adapted 
to the BBS model as a cyclic system. We start with the necessary first step, the 
solution of the associated auxiliary problem. In a second step we obtain the eigenvectors 
and eigenvalues of the periodic system by Baxter equations. The conditions which 
ensure that the Baxter equations have non-trivial solutions are formulated as truncated 
functional equations. Section 4 gives a description of local spin operators in terms of 
global elements of the monodromy matrix. Starting with Section 5 we restrict ourselves 
to the case N — 2, for which the BBS model becomes a generalized 5-parameter 
plaquette Ising model. In section 6 we further specialize to the homogenous case 
and then to the two-parameter Ising case. Periodic boundary condition eigenvectors 
are explicitly constructed. Section 7 is devoted to our main result, the proof of the 
factorized formula for Ising spin matrix elements between arbitrary finite-size states. 
This is shown to agree with the Bugrij-Lisovyy conjecture. In section 8 we give an 
analogous formula for the Ising quantum chain in a transverse field. Finally, section 9 
presents our Conclusions. Large part of this paper relies on our work in [25, 26, 27]. 
Sections 4, 2.2 and 6.2 give new material. 
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2.1. The inhomogenous BBS-model for general N 

We define the BBS-model as a quantum chain model. To each site k of the quantum 
chain we associate a cyclic L-operator [29, 30] acting in a two-dimensional auxiliary 

space 

/ 1 + AxfeVfc, Au^^(afe - 6feVjt) \ 
Lk{X)= I fc = l,2,...,n. (1) 

\ Ufe(cfe - dkVk), >^akCk + Vkhdk/ >ik J 

A is the spectral parameter, n the number of sites. There are five parameters x^, 
Qfe) bk, Cfe, dk per site, and are elements of an ultra local Weyl algebra, obeying 

UjUfc = UfcUj , VjVk = VfcV^- , UjVk = uj^^'''VkUj , uj = e^'''/^, = vf = 1. 

At each site k we define a A^-dimensional linear space (quantum space) Vk with the basis 
|7)fc, 7 G Zjv, the dual space with the basis jk(7|, 7 € Zjv, and the natural pairing 
fe(7'|7)fe — ^'r','y In ^fc and the Weyl elements Ufc and act by the formulas: 

Ufc|7)fc = cu^|7)fc , Vfc|7)fc = I7 + k{l\uk = k{l\ , k{lWk = k{l - M ■ (2) 

The monodromy Tn{X) and transfer matrix t„(A) for the n sites chain are defined as 

r„(A) = Li(A)---MA)=(^^j^J ^j^lj, t„(A)=trr„(A) = A„(A) + i?„(A). (3) 

This quantum chain is integrable since the L-operators (1) are intertwined by the twisted 
6- vertex i?-matrix at root of unity 

( \ -uv 0\ 

cj(A-t/) A(l-tu) 

1/(1 -a;) \-v 

\0 \-ujv J 

R{X,u) (A)Lf (^) = (^) (A) R{X,u), (5) 

where ^^^A) = Lk{X)®I, 4^^(A) = I®Lfe(A). Relation (5) leads to [t„(A), t„(/x)] = 0. 
So t„(A) is the generating function for the commuting set of non-local and non-hermitian 
Hamiltonians Hq, . . . , H„: 

t„(A) = Ho + HiA + • • • + H„_iA"-^ + H„A". (6) 

Prom (5) it also follows that the upper-right entry Bn{X) of T„(A) is the generating 
function for another commuting set of operators hi, . . . , h„: 

[Bn{X), BM] = 0, Bn{X) = hiA + h2A=^ +■■■ + h„Ar (7) 

Observe that Hq and H„ can be easily written explicitly in terms of the global Zjv-charge 
rotation operator V„ 



R{X,u) 



(4) 



n 



Ho = 1 + V„ JJ , H„ = JJ Ofe Cfe + V„ J]^ Xfe , V„ = V1V2 • • • v„. (8) 

k=l ''^^ k=l k=l 
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Here wc shall not explain the great interest in the BBS-model due to a second 
intertwining relation in the Wcyl-spacc indices found in [29] and the related fact that for 
particular parameters the Baxter Q-operator of the BBS-model is the transfer matrix 
of the integrable Chiral Potts model, see [29, 33, 38]. We will also not discuss the 
generalizations of the BBS model introduced by Baxter in [39], and not explain how (1) 
arises in cyclic representations of the quantum group Ug{sl2), see e.g. [23, 40, 41]. 

The transfer matrix (3) can be written equivalently as a product over face 
Boltzmann weights [28, 33]: 

tn(A) = Utl -1) 7jfc) Ik, Ik-i) with the face Boltzmann weights 

w^r(7fe-i, 7k, iLi, ik) = J:1,_,=o^'^'-'^'''-"'-'^ 

X iF^_i(7,_i - 7^_„ m,_i) i^'(7fe - 7^, m,^,) (9) 

where e {0,1} and F^(A7,mfe) = F^'(A7,mfe) = if A7 7^ {0,1}, and the non- 
vanishing values are 

""^^(x. -W-)' "^'^{l -dkJ-k)- 

The vanishing of F^(A7, nik) and F^'(A7, nik) for A7 7^ (0, 1} means that the vertically 
neighboring Zjv-spins cannot differ by more than one. The equivalence to the transfer 
matrix defined by (1) and (3) is seen writing the matrix elements of (1) as 

ilklLki^U.^mM = a;-'=-^^-'=^'=A^^-^— ^^^7^ - 7.,m,_0F^(7^ - 7.,^^). (11) 
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Figure 1. Illustration of the two versions: left we see the W of (9) indicated by full 
lines, whereas (3) the of (1) arise if we look at the lattice formed by the dashed 
lines in the left figure and the dashed rhombus shown at the right. The lattice is built 
by Zjv-spins on the full lines and Z2-spins in the centers. 



2.2. Homogenous BBS-model for N — 2 

The integrability of the BBS-model is valid also if the parameters x^, a^., . . . , vary 
from site to site and the construction of eigenvalues and eigenvectors can be performed 
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for this general case. However, in order to obtain compact explicit formulas for matrix 
elements, we shall often put all parameters equal: >Ck = >i, ■ ■ ■ , dk — d and call this the 
homogenous model. In [42] it has been shown that for N = 2 the general homogenous 
BBS-model can be rewritten as a generalized plaquette Ising model with Boltzmann 
weights 

Vl^(cri,(T2,(73,(74) = ao(l + I]i<j<j<4 Oij + 04(71(72(73(74^ , (12) 

subject to the free-fermion condition — 012034 — 013O24 + 014O23. 

For N = 2 the Weyl elements can be represented by Pauli matrices. Fixing h—\ 
the L-operator becomes 

\-r\al Xal{a-bal) 
al{c—dal) Xac+a^bd 

degenerating at A = b/a: 

1 + b/aa^ 



Lk{b/a) 



1, b 



The matrix elements of the corresponding transfer-matrix are 

n 

{{a'}\tn{b/a)\{a}) = Yl{5a,y^{l + bcak-ia'^) + S^^,_^'^b/a{l - adak-ia'^)) , 



k=l 



where {a} — {ai, . . . , (7„} and {a'} — {a[, . . . , a!^} are the values of the spin variables 
of two neighboring rows, ak — (— l)'"'*, a'^ — (—1)'''*^ £ {+1, —1}, and the identifications 



a'f. are used. 







< 




Ky 


/ Ky 








<^k 



Figure 2. Transfer-matrix for the triangular Ising lattice. Solid lines show the 
interaction between spins. 



The matrix elements of the transfer-matrix of the Ising model on the triangular 
lattice (see Fig. 2) are 



({cr'}|tA|{cr}) = W exp{K^ ak-iCTk + Ky akcr'^ + Kd ak-ia'i.] 



(13) 



k=l 
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The /c-th factor of this product, taken at — cf'^ is 
exp(Xj,) ex.Y>{{K^ + Ka) Ok-i cr^) = 

= ex^{Ky) cosh(ir^ + Ka) (1 + tanh(i^3; + Ka) (Tk-i cr'^) , 

and at ak = —cr'^ is 

exp{-Ky) exp{{Kd - Kr,) ak-i (j'k) = 

= exp(-A'j,) cosh(i^<i - K^) (1 + tanh(X<i - K^) ak-i cr'^) . 

Now it is easy to compare the transfer-matrices t„(6/a) and tA : 

tA = exp(ni^,) cosh"(i^, + K^) K{b/a) , exp{-2Ky) ''''tl^'' I = V« , 

cosn^Ad + Kx) 

tanh(ir2; + Ka) — be, tanh(ira; — Ka) — ad. 

Although we considered t„(A) at the special value of the spectral parameter X = b/a, 
the transfer-matrix eigenstates are independent of this choice of A. So the eigenstates 
of the transfer-matrix of the Ising model on the triangular lattice appear as eigenstates 
of the general homogenous BBS-model for A?" = 2 (the parameter x and one of the 
parameters a, d in the case of homogeneous periodic BBS model can be absorbed 
by a rescaling of the other parameters and using a diagonal similarity transformation of 
the L-operators) . The formulas for them will be given later. Unfortunately, factorized 
formulas for the matrix elements of the spin operator in this general case have not been 
found. There are only two special cases for which such formulas are available: 

• Row-to-row transfer-matrix for the Ising model on the square lattice: 

a = c, b^d: ATd = 0, e"^^" = 6/a , tanh = ab. (14) 

This case will be the main object of our attention. It is the most general case 
where we have the factorized formula for the spin operator matrix elements found 
by Bugrij and Lisovyy. 

• Diagonal-to-diagonal transfer-matrix for the Ising model on the square lattice: 

a = c, b^-d: K^-^O, e'^^^'^b/a, tanhKa^ab. (15) 

It is known [43] that such transfer-matrices with different parameters Ky — L, 
Kd — K (and corresponding a, b) constitute a commuting set of matrices having 
common eigenvectors provided 

smh2K sinh2L = — = -j (16) 

1 — 0^ k 

is fixed. Thus in this case the eigenvectors depend on k' only. Therefore, in order 
to find the eigenvectors and the corresponding matrix elements of the spin operator 
it is sufficient to fix a = c = l/(k')^/2 ^j^d b = d = and so to obtain a special 
case of the formulas for the row-to-row transfer-matrix of the Ising model on the 
square lattice. Note that we get [27] the same matrix elements in the case of the 
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quantum Ising chain in a transverse field with strength k' because the corresponding 
Hamihonian commutes with the transfer-matrices having the same k'. Another 
remark: with the restriction a = c, b = —d, x = 1, the transfer-matrices commute 
among themselves at independent values of two spectral parameters: A and the 
parameter which uniformizes (16) (a parameter on the elliptic curve with modulus 
k', see [43]). 

3. Separation of Variables for the cyclic BBS-model 

3.1. Solving the auxiliary system (7): Eigenvalues and eigenvectors of Bn{X). 

We start giving a summary of the SoV method as applied to the general inhomogenous 
Ziv-BBS model [25]. The aim is to find the eigenvalues and eigenstates of the n-site 
periodic transfer matrix t„(A) of (3), and the idea [19, 20, 21, 22] is first to construct 

a basis of the A^"-dimensional eigenspace from eigenstates of i?„(A), see (7). This can 
be done by a recurrent procedure. Then the eigenstates of t„(A) arc written as linear 
combinations of the A) -eigenstates. The multi-variable coefficients arc determined 
by Baxter T — Q- equations which by SoV separate into a set of single- variable equations. 

Prom (7) the eigenvalues of Bn{\) are polynomials in the spectral variable A. 
Factorizing this polynomial, for n > 2 we get 

n-l 

S„(A) = AAo n (-^ - Afe) I*a); a = {Ao, Ai, . . . , A„_i}, (17) 

fe=i 

where Ai, A2, . . . , A„_i arc the n — 1 zeros of the eigenvalue polynomial and Aq is a 
normalizing factor. We can label the eigenvectors by A. 

An overview of the space of eigenstates of -B„(A) is easily obtained using the 
intertwining relations (5). It follows from (5) that the operators An{X) and -Dn(A) 
of the monodromy (3), taken at a zero A = Aj^, are cyclic ladder operators with respect 
to the kth component of A in |^a)- To see this consider e.g. the intertwining relation 

(A - u;fj,)AMBn{fJ') = u;{X - ij,)BMM^) + - u;)AMBn{X) W 
which is a component of (5). Fixing A = A^, A; = 1, . . . , n — 1 , in (18) and acting on 
"ifx, the last term in (18) vanishes and we obtain 

5„(/.)(A„(Afc) I^a)) =/iAo(/x-a;-%)n3^fc(/^-A.) (A„(Afc) |^a)).(19) 
This means that 

An{Xk) = <Pk ■ |*Ao,...,a,-iA„...,A„_i> ■ (20) 

Later we shall give an explicit expression for the proportionality factor (f^. Similarly, 
from another component of (5) and with another factor (p^ we get 

Dn{Xk)\-^x) = ^fc -1*0,-^0, ...,a,A,,...,A„_i>- (21) 

Furthermore, acting by (18) on |^a) and extracting the coefficient of A"+^//" we get 

Vnl^'A) = |*u,-iAo,Ai,...,A„_,)- (22) 
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Assuming generic parameters in such that all proportionality factors are non- 
vanishing, by repeated application of An{\k), Dn{\k) and V„ to any eigenstate we 
span the whole A^"-dimensional space of states. Later, when we give explicit expressions 
for ipk and (fk we can check whether these factors can vanish. 

So, if for a given set of parameters ak,bk, Ck, dk, Xk, {k = 1, ...,n) there is an 
eigenvector with the eigenvalue polynomial determined by the zeros A, then there are 
also eigenvectors to all eigenvalue polynomials determined by the zeros 

{Ao ... , A„_i uj""--'} with p„ = (p„,o, • • • , Pn,n-i) e (Z^)". (23) 

Let us therefore write the zeros as 

K,k = -rn,feu;^"'N (24) 

where for n fixed, the n real numbers r„ ^ are determined by the 5n parameters ai, . . . , x;. 
For fixed parameters the iV" in all following calculations we shall label the eigenvectors 
by the p„ instead of our previous Xn,k- For given parameters, the set of the eigenvalues 
is determined by the Tn^, . . . , The eigenvalue equation for Bn{X) becomes 

n-l 
k=l 

In order to calculate the rn,k in terms of the parameters, we don't need the full 
quantum transfer matrix and the L^-operators involving the Weyl variables. Rather, by 
the following averaging procedure [23] 

0{X^) = (0(A^))= Usez.Oiu^^). (26) 

we associate to a spectral parameter dependent quantum operator 0{X) a classical 
counterpart O(A^) . We define the classical BSS model by the L-operator >C^(A^) 

[Loo) {Lo,)\ (I- e.^A^ -6A^(a^ - h^) 




(27) 



1) . Analogously, we define the classical monodromy Tn by 



r. = A(A")..(A'',....„(A'')^(^;;j; ^;f;;) (28, 

Proposition 1.5 of [23] tells us that the classical polynomials An{X^), Bn{\^), Cn{X^) 
and Vn{\^) are the averages of their counterparts in (3): ^n(A^) = {An{X)), etc. 
So for n > 2 we have 

m— 1 

B^(A^) = (-e)-A^<o n (^"^ - '<s)- (29) 

It is easy to derive [25] a three-term recursion which expresses Bm{X^) in terms of 
B„,-i{X^) and Brn-2{X^). Using the initial values Bi{X^) = -eA^rf; Bo{X^) = 
and defining rf^ = af' — , this gives a (n — l)t/i-degree algebraic relation for the r^^. 

For the homogenous model (the constants are taken to be site-independent) this 
can be replaced by just a quadratic equation, see the Appendix of [25]. 
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3.2. Solving the auxiliary system: Explicit construction of the eigenvectors of Bn{\). 

The stepwise construction of the eigenvectors, starting with one-site, then two-site as 
hnear combination of products of two one-site eigenvectors etc. is tedious because we 
have to go to 4 sites before the general rule emerges. 

Let us start finding the one-site right eigenvectors \il^p)i of -Bi(A) as linear 
combination of spin states 1 7)1, 76 Zjv, writing 

= X] ~ P^^N- (30) 

Applying on the left Bi from (1) and on the right (25), we demand 

Aur'(ai-6ivi) J2 - P) l7)i = An.ou;-"^-" Yl " (^1) 

7 6 Zat 7 e Zjv 

Applying (2) and shifting the left hand summation for the term with I7 -|- we get 

(oi - ri^oi^'^''') Wp{-f - p) = biWp{-f - p -1). (32) 
This is a difference equation for the function ^^(7) [44]: 

^7^ = r^^; ^p(0) = l; TeZiv, (33) 

Wp['y — 1) 1 — ui^ X 

where we have put y = b-i/a-i, ri o = xai and chose the initial value Wp{0) = 1. The 
cyclic property ^^(7) = Wp{j + N) imposes the Fermat condition x'^ + y'^ — 1 on 
the two-component vector p — {x, y). We indicate p as a subscript on the functions 
Wp{'y). We shall consider the case of "generic parameters", so in particular we exclude 
the case — b^ = 0, and the "superintegrable" case 

Qk = ou~'^bk = Cfe = dfe = Xfc = 1, (34) 

since in the latter cases degenerations appear. 
We write the analogous left eigenvector as 

^(^pl= E ^^(7-p-l) ^^^^ ^^^^ 

with the same functions Wp{'y), just now p = {ri^/ai, u}~^bi/ai). The Fermat vector 
dependent functions ^^(7) play an important role for cyclic models. They are root-of- 
unity analogues of the g-gamma function. 

By a similar calculation, the two-site eigenvectors are found to be: 

E^-(P2,0+P2,l-Pl)(p2,0-P2) 
-( T\ 7 7T IV'pi)! ® |V'P2)2- (36) 

where p2,o = (3^2,0,2/2,0), P2 = (^2,^2) and 

ri r2 . r2 _ 62(^2 n .„„s 
X2,o = 02 C2 , 2/2,0 = >ii , X2 = , 2/2 = ■ (37) 

r2,o r2,o r2,or2,i ^2 ^^2,0^2,1 

The condition that p2 and p2 are Fermat vectors determines r2o and r2 i- 
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The explicit formula for both the left- and right eigenvectors of Bn{\) for general 
number of sites n has been proved by lengthy induction and is given in [25]. A by- 
product of these calculations are the formulas for (p^ and (fk introduced in (20), (21): 

n-2 



s=l 



n-1 



Dn{Xn,k)\'^p„) = Mp'u) l*pr--^)' 'fiM = " r.n-2 nk H ^rn{Xn,k) , (39) 

'"-1 lls=l yn-l,s m=l 

Fn{X) = {bn +UJan>i:nX) (Ac„ +dn/Xn)- (40) 

On the left of (38) and (39) the eigenvectors of -B„(A) are labeled by the vector 

Pn = {pnfi, Pn,n-l) G (^Jv)"- (41) 

p^'^ denotes the vector p„ in which p^^^ is replaced by pn,k ± 1: 

Pn'' = (Pn,0,---,Pn,fc±l,---,Pn,n-l), /c = 0, 1, . . . , 71 - 1, (42) 

Tn = rnflr„,i . . . r„,n_i and Pn = J2k=oPn,k- (43) 
p'^ denotes the vector p^ without the component p„^o- 

= (Pn,l,---,Pn,n-l)e(Ziv)"-^ (44) 

The yn-i,s components of a Fermat vector Pn'\,s = (^n-i,s' yn-i.*) <^cfincd by 
^n^i s = rn,k/rn-i,s, sce Section 2.4 of [25]. The ^^(A) which appears in (38) and (39) 
is a factor of the quantum determinant: 

n 

An{uX)Dn{X) - Cn{ijX)Bn{X) = V, • J] F^{X) , (45) 

m=l 

Prom (1) we can read off directly the A°- and A"-coefficients of the polynomial A„(A): 

A^{X) = 1 + ... +xiX2...x,VA". (46) 

Then using (38), the general action of An{X) on Bn eigenvectors can be written as an 
interpolation polynomial 

n— 1 / , \ n—1 

^n(A)i*p„) = n(i-i— ) i*pj+Axi---x„ n(^-^n,^)i*pr)+ 

s=l ^ '^"'^^ s=l 

Considerable effort is needed to present the norm of an arbitrary state vector 
|^p„) in factorized form, since multiple sums over the intermediate indices have to 
be performed. The norms are independent of the phase Pn,o and their dependence on 

Pn is: 

m 1^ \ = — = — (AS) 
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The normalizing factor C„ is independent of p„ and can be written recursively [26]. The 
two lowest values are: 

<^ \yij ^ \y2y2y2,oJ 

3. 3. Periodic model: Baxter equation and truncated functional equations 

In the auxiliary problem we looked for eigenfunctions of -B„ does not commute with 
V„ (8), see (22): V„|^p,J = |^p+o). Now we are looking for eigenfunctions of t„ which 
commutes with V„. By Fourier transformation in p„ o we build a basis diagonal in V, 
where the Fourier transformed variable p E is the total ZAr-charge.: 

I^px) = Ep„,oez.^-''-'""°l*pJ, V„ = o;^ \^,,^J . (50) 

We now write the eigenfunctions |$p,e) of t„(A) as linear combination of the |^p,p;). 
The eigenvalues of t„(A) on these states are again order n polynomials in A: 

t„(A)|$,,E) = (£;o + £;iA + --- + K-iA"-^ + KA")|$p,E> (51) 
Since the values of Eq and En can be read off immediately from (8): 

^0=1+ ^''nm=l Kid„i/ X„„ En = YIZ,=1 (^rnCm + IlrUl ^rn , (52) 

wc combine the remaining coefficients into a vector E = {£^1, . . . , En-i} , and label 
the eigenvectors just by the charge p and E: 

t„(A) = tn{X\p,E) |$,,e), |*p,e) = E ^''(PnlAE) |^,,p,). (53) 

P'n 

Now, in order to achieve SoV of the multi-variable functions Q^, we split off from 
Q^{p'n\ p, E) Sklyanin's separating factor: 

Q^(p;i P, E) = n^^i (^"-^) . (54) 

li.,.'=i W „y[Pn,s - Pn,s') 

We shall not give the detailed calculation and just indicate the main mechanism. We 
express t„(A) as an interpolation polynomial through the zeros A„^fe of S„(A): 

(^„(A) + D„(A))|«,,^.> = n (1 - 1^) + a -B„n(A - A„,,)| + 

+ E (U^^)!^A^^M\K^t') + -'M)\^,.^,^))- (56) 

When we evaluate (55) successively at the n — 1 values A = \n,k, k = 1, . . . ,n — 1, the 
terms on the right of the first line of (55) do not contribute. Due to the Sklyanin-factor 
the brackets involving the differences Xn,k — Xn,s are made to cancel, leading to SoV. 
This results in the n — 1 single- variable Xn,k Baxter equations (A; = 1, . . . , n — 1) 

tn{Xn,k\P^ E) Qf{Pn,k) = A+(A„,fc) Qf{Pn,k + 1) + A^(wA„,fc) QUpn,k " 1) ■ (56) 
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Starting from the left eigenvectors the analogous left Baxter equations are 

tn{\n,M E) Q\{Pn,k) = o;""' A" (A„,fe) Q\{pn,k + 1) + o;'"" A+(u;A„,fe) Q\{pn,k - 1) ,(57) 

where we abbreviated 

n-l 

A+(A) = (u;7x;t) (A/o;)^-" J] ^-(A/o;) , A-(A) = Xfc (A/u;)"-' i^n(A/a;) ■ (58) 

m=l 

Xfc collects several factors (partly arising from ip^ and ipi^) determined by constants 
Xfc, Ofc, . . . , dfc alone. Now note that the left hand side of (56) more explicitly reads 

where the E are unknown and have to be determined from the system of homogenous 
equations (56) together with the n — 1 functions Q^{pn,k)- In order to have a non- 
trivial solution, the coefficient determinants have to be degenerate. Fix a k, then 
from the determinant we may get one relation among Eq, . . . , En. All n — 1 systems 
for different k should be sufficient to determine all components of E. Fortunately, 
the condition for non-trivial solutions to (56) can be written as well-known truncated 
functional equations: 

Define t(^)(A) = ^(A)t and construct a fusion hierarchy [45, 33] by setting 
t(o)(A) = 0, T«(A) = 1, and 

r(^-+i)(A) = T(')(a;^-^A)r(^)(A) - ^''^(c^'"'A)r(^-^)(A), j = 2, 3, . . . , A^, (59) 

ziX) = a;-''A+(A)A-(A) = ni=iFm{XM. (60) 

Then it can be shown [25] that if r*^-^"^^) (A) satisfies the truncation identity 

^(iv+i)(A) - uj''z(X)t^''-'\ujX) = A(A^) + P.(A^) (61) 

with ^„(A^) + Vn{X^) given in (28), then the system (56) has a non-trivial solution 
for all k. This truncated hierarchy can be used to find the transfer matrix eigenvalues 
[46, 47]. In our construction we have even more: for every solution of (59), (61) we can 
construct an eigenvector. 

4. Action of and on eigenstates of -B„(A) 

Our main aim is to calculate matrix elements of the local operators uj. and between 
eigenstates |$p,e) of t„(A) . Since we know how to get these states from the -Bn(A) 
eigenstates (53), (54), (56), we first set out to find the action of the local operators on 
the |^E'p„). Since we built our auxiliary states successively from one-site to n-site, the 
formulas will not be symmetric between e.g. Uj and with j ^ k. 
For u„ we can calculate its action directly. Starting from 

I This definition in [29] is the origin of calling the BBS model the r*^^-* -model 
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we get the formula for the action of u„ on one-site eigenvectors: 



to- 



Pn 



Unl^pjn = K|V'p„)n - &n|V'p„+l)n) • (62) 



Using then the exphcit recursion formula relating |^p„) to |^p„_i) one finds [26]: 

UT, \ l,T, \ ^n>^l>^2 • • -^^n-l \ , , ao\ 

*Pn) = ~ ^Pn ZT^^ ^p+o + 63 

We shall derive this result in a simpler way expressing the local operators and in 
terms of the global entries and i?„ of monodromy matrix, taken at particular values 
of A. There is a well-known method elaborated by the Lyon group [48]. However, this 
method requires the fulfillment of the condition i?(0) = P with R the quantum i?-matrix 
intertwining two L-opcrators in quantum spaces and P the permutation operator. This 
requirement is fulfilled for the cyclic L-operators only at special values of parameters 
where the i?-matrix is the product of four weights of the Chiral Potts model [29]. 
Another requirement regards the possibility to obtain such a it!-matrix by fusion in 
the auxiliary space of the initial L-operator. This requirement can not be fulfilled for 
the cyclic L-operators (1) because the fusion in the auxihary space [49] gives L-operators 
with the highest weight evaluation representations of the corresponding quantum affine 
algebra, but we need cyclic type representation in the auxiliary space. 

We will use an idea borrowed from a paper of Kuznctsov on SoV for classical systems 
[50]. What we can do is the following: Consider the inverse of the operator Lfc(A): 

/ a; A Ofc Cfe Vfc 6fe 4/xfe -Au^^^ (ofe - 6fe v^) \ 

L-^\\)= ■ (det,L,(A))-S (64) 

y -u;Ufe(cfc - 4vfe), 1 +a;AxfcVfe ] 

where 

det5Lfc(A) = VfeFfc(A) , Ffe(A) = (6^ u\a},K^{\c^ 4/>^fe) ■ 

The expression for Lj^''^(A) is singular at zeros A^ = —hk/{uja};>Ck) and A'^ = —dk/{ck>^k) 
of Ffe(A). Of course, 

T„_i(A) = T^{\)L-\X). (65) 

Therefore at the zeros of L'n(A) the left-hand side is regular in A and the right-hand side 
also has to be regular. At A = A^ = —hn/ (coan^Cn) we get 

Hence we have a formula for u„: 

u„ = A;a„S-^(A;)^„(A;). (66) 

Prom the condition of the regularity of the right-hand side of (65) at A = A^ = 
-4/(c„x„) we get 

An{K){-K)K\(^n-bny^n)+B4X':^){l-dn/{u;Cn)Vr,) =0. 
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Excluding u„ by means of (66), we obtain the formula for v„: 

X (a,(a;) i?n(A::)/A:; - a^{\':,) bm.)/K)- (er) 

Using the RTT-relations following from (5), we can permute An and in (66) to get 
an equivalent formula 

u„ = a; a„ An{u;X'J B-\cjX'J . (68) 

Using (47) and (25) we get (63). 

Wc can also get the formulas for u„_i and v„_i. We express L^^(A) in terms of 
An{X) and -Bn(A) using (66) and (67). Now the formula (65) allows to find expressions 
for An-i{X) and i?„_i(A) in terms of A„(A) and Bn{\). Finally we substitute these 
expressions to (66) and (67) in which the indices n arc replaced by n — 1. This gives us 
expressions for u„_i and v„_i in terms of A„(A) and -B„(A). The described procedure 
can be iterated to express the local operators and Wk in terms of An{\) and -B„(A). 
For example, the result for u„_i is: 

U„_i = UJ>^^_^an-l(Ar,{uj}^^_^){u'^M^_^anCn + ^nKdn/ Xn) - 5n(t^A^_i)u;U„(c„ - rf„V„)^ 

X (^-An{uj>!^_^)u>!n_in~^{an - bnVn) + -B„(u;A^_i) (1 + cu^A^.iX^Vn) j , 

where u; A^_]^ = — 6„_i/(a„„i x„_i) and the expressions (68) and (67) for u„ and v„ 
have to be substituted. It gives the action of u„_i on |^p„). We see that the formula 
gets quite involved. However, ui can be easily expressed in terms of and B^. 

u, = i D„ (-^) B-^ (-^) . (69) 

For our purpose of finding matrix elements of spin operator between eigenstates |$p,e) 
of homogeneous t„(A) we can choose any spin operator because they all are related 
by the action of translation operator having the same eigenstates |$p,e)- In what follows 
we consider matrix elements of the spin operator u„ because the corresponding formula 
for the action (63) is the simplest. 

At the end of this section we would like to mention some similarity of our formulas 
with the formulas from the paper [51], where the local operators for the quantum Toda 
chain are expressed in terms of quantum separated variables with the use of a recursive 
construction of the eigenvectors [22] . 

5. The general inhomogenous N — 2 BBS-model 

In the N = 2 case we have two charge sectors p = 0, 1. Following the language of 
e.g. [14, 17, 24] the sector p = will be called the Neveu-Schwarz (NS)-sector, and 
p = 1 the Ramond (R)-sector. We are going to show that the spin matrix elements 
can be written in a fairly compact, although not yet factorized form (85), (86). The full 
factorization will be achieved later for the homogenous Ising case. 
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5.1. Solving the Baxter equations and norm of states 

Let us fix an eigenvalue polynomial t(A|p, E) of t(A) corresponding to a right eigenvector 
|^p,e) (since in the following our chain will have the fixed length n we often shall skip 
the index n. Also sometimes we shall suppress the arguments p, E in i). 

In order to find |$p,e) explicitly we have to solve the associated n — 1 systems 
(A; = 1, 2, . . . , n — 1) of (right) Baxter equations: 

t{-rn,k) g^(0) = (A+(-r„,fe) + A^(r„,fe)) Q^(l), 

t(r„,fc) QUI) = (A+(r„,,) + A,(-r„,,))g^(0). (70) 

Since t(A|p, E) is eigenvalue polynomial, the functional relation (61) ensures the 
existence of non-trivial solutions to (70) with respect to the unknown variables Qfc (0) 
and Qfc (1) for every k — 1, 2, . . . , n — 1. In the N — 2 case, this means that for every 
k we have one independent linear equation (in case of degenerate eigenvalues, possibly 
no equation). In the case of generic parameters, both hand sides of each equation will 
be non-zero. So, fixing Qf{0) — 1 we obtain two equivalent expressions for Qf{l): 

^Rf. N ^ t{-rr,,k) ^ A+(rn,fc) + A-(-rn,fc) , . 

^ A+(-r„,,)+A,-(r„,,) t(r„,,) " ^ ^ 

Analogously from the left-Baxter equations, fixing Q^(0) = 1 we obtain 

^L^.x ^ {-ir-%-rn,k) ^ A+(-r.,,) + A-(r.,,) 
^^^^ A+(r„,,)+A-(-r„,fc) (-l)"-4(r„,,) " 

Since for generic parameters t(r„ fc|p, E) ^ these explicit formulas give 

QkiPn,k) QUPn,k) = (-l)''».'=("-i) ^((-1)''-'= r,,fc)/t(r,,,) . 

To get the periodic state, we have to insert the Skylanin-separation factor (54). Now 
for N = 2 the functions Wp are simple: 

^,(o) = i, ^p(i) = YT^ = ^> (^^(1))' = irf • ^^^^ 

In the Sklyanin factor we have to use the Fermat point p^'™ = (a;^'^, 2/^'P) defined by 
the coordinate x^'J" = rn,m/ fn,i ■ Here it can be expressed it in terms of x^'J" only and 
we get 

('^P.eI'^P.e) ^ nroi(-l)^"''^^-'"(^n,^ + rn,l? Hill Q\{Pn,k)Q't{Pn,k) .73. 

We can normalize to a convenient reference state. For the moment, simple formulas arise 
if for the normalization we chose the auxiliary state |^o,o) where = (0, 0, . . . , 0). From 
(48) we get 



' n,m + 



(74) 



Ising model spin matrix elements from Separation of Variables 16 

Combining all these formulas we get for the left-right overlap of the transfer matrix 
eigenvectors of the periodic BBS model a,t N — 2: 

(^P.E I '^p,e) ^ U':<rLirn,m + rn,l) ^ 111=1 (-^Y"'' ^((-l)^"-'r„,0 

This formula is not yet very useful since from (53) it contains the summation over 
the n — 1 Z2-variables p'^ defined in (44). However, in [26] it is shown how to perform 
this sum explicitly, and the fully factorized result is 

(^p,E I ^p,e) _ r,n-l~> IYl<mi'^n,m + ^n,;) yr , . . 

(*0,o|*0,o) Wk=l\Al=l irn,l + l^k) 

where — /^.j are the zeros of the eigenvalue polynomial of t(A|p, E): 

n 

t(A|p,E)|$,,E) = A n + A^O I^p,e). (77) 

i=l 

We don't specify the factor A, since in the following it will cancel. 

5.2. Matrix elements between eigenvectors of the periodic N = 2 BBS model 

In (63) we obtained the action of u„ on an eigenvector |\E'p„) of Bn{\): the result is 
a linear combination of the original vector plus a sum of vectors which each have one 
component of p„ shifted. In order to get the matrix elements of u„ in the periodic 
model, using (50) we first pass to charge eigenstates (^p^pj^l, \^ p^pij): 

(*P,P'J = (*oxl + (-)''(*i,P'J , I^P-Pi.) = l*o,p'J + (-)1*i,p'J- (78) 

Since uj — —1, u„ anti-commutes with V„ so that only matrix elements of u„ between 
states of different charge p can be nonzero. In the following we shall chose the right 
eigenvector from p = 1, then the left eigenvector must have p — Q (the opposite choice 
gives a different sign in (79)). Using (63), we find 

(^oxl^ox) ^-.0 



(80) 

Of physical interest are the matrix elements between periodic eigenstates. To get these 
we have to form linear combinations determined by the solutions of the Baxter equations: 
Recall (53): |$p,e) = 2^(PrJP)E) and the corresponding left equations. 

Let ($o| be a left eigenvector of the transfer- matrix t„(A) with p = and |$i) be 
a right eigenvector with p = 1 (often suppressing the subscripts E, E'): 

($o,E'|t(A|0,E') = i(°nA) (*'o,eI , t(A|l,E)|$i,E) = t^'\\)\^i,^). (81) 
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Let Q^^^\pn,k) and Q^^^\pn,k) be the solutions of Baxter equation corresponding to 
these two eigenvectors. After some simphfication wc get for the matrix elements (keeping 
the normalization by the auxiliary "reference" state): 

Y.N(p-) {r,(p') (^(-1)'' - "'""•^^"-'''" ) + £ , (82) 



(^0 


\a 


■z 1 
n 




(^o,o| 


*0,0 > 



where 



n— 1 n— 1 



, n— 1 

X (l-^)^f^. (84) 

with Tk = r„,fc, pfe = A; 0, 1, . . . , n - 1, 

i^fe = -rfc(-l)'''= , f = ro ri • • • r„_i and p' = YTk^i Pk- 

The origin of the different terms in (82) is: the sum over p' comes from (53), M{p') is 
the normalization factor from (74). The terms at Ro{p') arise from the first hne of (63): 
the shift in p„ o affects the charge sector only. The sum over k and expression for Rk{p') 
come from the second line in (63). Now, the sum over k can be performed. Indeed, as 
shown in [27], using the Baxter equations, some cancellations take place and (82) can 
be written as 

= ^ E mp')RMmp') (85) 

( ^co I *o.o ) 2r„^^^,_, 

with 

^"^^ " nr/(-c.+(-i)-n) + nL-i(c+(-i)«n)' ^" ^''^ 

Despite the simple appearance, for the general inhomogenous N = 2 BBS-model, 
performing the sums over the Z2 variables explicitly seems to be a presently hopeless 
task. However, for the homogenous Ising model we shall show this to be possible. 

6. Homogeneous N — 2 BBS-model 

6.1. Spectra and zeros of the Bn- and tn-eigenvalue polynomials 
We now speciahze to A?" = 2 taking all parameters site-independent ( "homogenous" ) : 
Om^o, hm^h, Cm^c, dm ^ d, = X, = r, >C^(A^) = >C(A^), Vm. (87) 
Then the classical monodromy is 



(£(A^))". (88) 



Ising model spin matrix elements from Separation of Variables 18 



Consider trace, determinant and eigenvalues x± of C: 

5{\^) = det £(A') = (fcVx" - A'a') {dF - X^cW) = F{\) F{-X), (90) 



t(X^) = tr£(A') = 1 + ^ - X\k^ + a^c^), (89) 



a;± = i(r ± - 4 5), F(A) = (6 - axA)(Ac + rf/x). (91) 

From the matrix >C(A^) we obtain 

B^(A^) = -A^ {a' - b') {xl - xl)/{x^ - x_), (92) 
so that the zeros of Bm are at x+/x- = e'"^"^"-^ with 

(f)n,s = 27vs/n, 6- = 1, 2, . . . , n - 1, 0. (93) 

Using = 4(5 cos^(0/2) and (89), (90) we can translate the zeros labeled by 4>n,s by 
a quadratic equation in A^ into zeros A„ 5. 

Now we solve the functional relations (59), (61) for the transfer matrix spectrum. 
Using (59) for j = 2 and eliminating r^^^ by (61) we get the functional relation 

t{X) t{-X) = {-inz{X) + z{-X)) + An{X') + V^{X') (94) 

which we shall use to find t{X). In terms of (89) and (90) this reads 

t{X) t{-X) = (-1)'' {51 + 5!!) + xl + xl. (95) 

where 5± — {b ± axA) {d ^ cxA); = ^{^^) — x+x^. Introducing q taking 

the n values 7r(2s + 1 — p)/n, s = 0, . . . , n — 1 , we can write (94) as 

t{X)ti-X) = {-irllie%+-TiX')+e-'^5.) = i-irll{A{q)X'- Ciq) + 2iB{q)X) 
q q 

^1^^ A(q) = a^c^ -2xaccosq + x^; B{q) = {ad - be) sin q ; 

C{q) = 1 - 2{bd/x) cosq + b'^d^K^ . (96) 

Factorizing the polynomial in A we get 

t{X) t{-X) = {-irllA{q) (A - s,) (A + s_,) (97) 

q 

«q = - D{q) = A(q) C(q) - S(q)^ (98) 



(fixing the sign of ^jD{q) requires a special convention, see [25]) and after some 
arguments we find the spectrum 

t{X) = (a-c'^ + (-l)''x'^) n,(A±«q), (99) 

where the signs are not yet fixed. Comparing the A-independent term in (51) 

t{X) = 1 + (-1)^ 6" d"/x" + £;iA + • • • + K-iA"~^ + A"(a"c" + (-l)''x"). (100) 

with the corresponding term in (99) shows that the number of minus signs in (99) must 
be even (odd) for the NS-sector p = (R-sector p — 1). 
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It is useful to introduce the following notion: The eigenvalue (99) with all +-signs 
is called to possess "no quasi-particle" excitations. Each factor labeled by q having a 
minus sign is said to contribute the "excitation of the q-quasi-momentum" . We shall 
accordingly label the minus signs by a set of variables o"q G Z2, where for unexcited 
(excited) levels q we put cxq = (cxq =1). So instead of (99), we shall write more 
precisely 

t(p)(A) = (a"c" + (-l)''x") (A + (-1)"^ Sq). (101) 

The corresponding eigenvectors have been considered for the inhomogenous case in 
Subsection 5.2. 



6.2. Functional relation for the diagonal-to-diagonal Ising model transfer-matrix 

In this subsection we specialize the results of the previous subsection to the case of 
the diagonal-to-diagonal transfer-matrix of the Ising model on a square lattice (15). 
So, we set a = c, 6 = — x = 1 and A = h/a. Let us calculate the ingredients of 
the functional relation (94). We have -Fm(A) = —{b — aX)^. Therefore due to (60), 
^(A) = (-1)"(6 + aA)2", z{b/a) = (-1)" (26)^", z{-b/a) = and the averaged 
L-operator (27) at A^ — b'^/ becomes 

Hence 

A(&V«') + ^^n(&V«') = trr„(6V«') = (l-&Va2)"(l-a262)n. 
Substituting these expressions into (94), we get the following functional relation 

t{b/a)t{-b/a) = (-l)''+"(26)2"+ (l-67a2)"(l-a262)". 
We want to compare this with the functional relation equation (7.5.5) in [43]: 

V(K,L) V{L + m/2,-K)C = (2isinh2L)" 7 + (-2isinh2X)" i?, 

where C is the operator of translation, R is the operator of spin flip V„ and 
V{K, L) is the transfer-matrix (13) with ^ 0^ Ky ^ L, = K, e~'^^ = b/a, 
tanhX — ab. Therefore V{K,L) — exp(nL) cosh^X t„(6/a). Similar analysis gives 
V{L -\- i7r/2, —K)C = i"exp(nL) cosh^X t„(— 6/a). Now taking into account that the 
eigenvalues of R are (—1)'' and 

, 4a6 , ^ a^ — b"^ exp(— 2L) , ,,,, , , 

2sinh2/s: = —, 2sinh2L = — , — — ^ = (l-a^b^)b a, 

l-a262' ab ' cosh^X ^ ^ ' 

we see that both functional relations are identical. 
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6.3. Ising model: Spectra and zeros of the S„(A)- and tn{X)- eigenvalue polynomials 
We now specialize further to the Ising case (14) as aheady advertised in Subsection 2.2: 

ttj = Cj = a, bj = dj = b, Kj = 1; Vj. (102) 

In the Ising case (102) the 2" eigenvalues of (101) with (98) can be written (2"~^ in each 
sector p = 0, 1): 



where the quasi-momentum q in each sector takes n values: 

2tx 

q = — m, m integer for p = 1 (R); m half-integer for p = (NS). (104) 

77/ 

Recall that we found from (100) that in the NS (R) sector, the eigenstates of t(A) have 
an even (odd) number of excitations: (^1)'^'' = (~1)^- 

For q = (this occurs for R-sector only) and q = tt we define 
b^-\ _ 62 + 1 

q = TT is in the R sector for n even. However, for n odd it is in the NS sector. The 
different presence of factors (A ±so) and (A isyr) in (103) for n even or odd often makes 
it necessary to consider the cases n-even and n-odd separately. In the following we shall 
reserve the notation Aq for Aq = (— l)°'''Sq and otherwise use Sq as defined in (103). 
The zeros A„^fe of the -B„(A) eigenvalue polynomial are determined by (93), (89), (90): 

r(A^fc) = 4cos2gn,fcF(A„,fc)F(-A„,fc), g„,fc=7rA;/n, A; = 1, n- 1.(106) 

Since now o n a a ^ 

F(A) = F(-A) = 6' - a' A'; t(A') = 1 + 6^ - (1 + a^) A^ (107) 

we get 



«o = -2 7 , = -^VT ■ (105) 



rn,k = \j{b^ - 2 b"" cos qn,k + l)/{a^ - 2 cos qn,k + 1) = s^^ , . (108) 
Observe that Sq and rn^k ^lay coincide. 

6.4. Ising model state vectors from Baxter equations 

In order to obtain the eigenvectors of t(A), we have to solve Baxter's equations. For 
our restricted parameters (102) we have F{\) = F{—X) and the left and right Baxter 
equations (57), (56) become identical. Omitting the superscripts L and R on and 
recalling An,fc = -(-l)''"''=r„,fc, pn,k = 0, 1 we obtain: 

^ ' "'"^ + {-KkT~^ Xk F{Xn,k) Qk{Pn,k + 1). (109) 

['^n,k) Xk J 

From (109) we get the following compatibility condition: 



t{-r^Mrn.k) = (-1)"-^ ( ^ (4,)n-i ^r + (-rn,.)"-^X.F(r„,.) ) , 
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if t{X) is an eigenvalue from the sector p. If (—1)'^ = (—1)''+^ then the quasi-momentum 
q = qn,k belongs to the sector p and for rn,k = Sq^ ,^ we have t{-rn,k)t{rn,k) = 0. This 
implies a relation not depending on a particular t{\) and its p: 

Xl rlt'^ - {-ir'^'F--\r^,,) . (110) 

Although the eigenvalue polynomial t{X) is known from (103), to solve (109) for the 
Qk{Pn,k) can meet a difficulty if tn{\n,k) vanishes or if, due to (110), the big bracket on 
the right of (109) vanishes. All this can happen and we have to distinguish four cases 
(we suppress n and write just = Tn,k and pk — Pn,k)' 

(i) (-1)'' = (-l)*": This is the easy case, since from (104) and (106) 
t'^i^k) 7^ and tP{—rk) 7^ 0, and we may normalize and solve 

Qk Qk (1)- 2xkrl-' F{r,) ■ 

The other three cases occur for (—1)^ = (—1)*^""^ : 

(ii) tP{rk) 7^ 0, tP{—rk) = 0: t^(A) contains a factor (A + r/j)^ (both q = ±qk not 
excited), we may normalize 

g^>) = i, g^'^(i) = 0. 

(iii) tP{rk) — 0, tP{—rk) 7^ 0: t^(A) contains a factor (A — Vkf (both q = ±5^ are 
excited), we cannot choose (^^'^(O) = 1, but we may normalize 

g^'^o) = 0, Q]:\i) = 1. 

(iv) tf{rk) — tf{—rk) — 0: ^''(A) contains (A^ — r^) (either q = +qk or q = — g^. is 
excited): A L'Hopital procedure, using a shght perturbation of (102) as described in 
[26], is required (to obtain eigenvectors of translation operator), leading to 

/ -|^yr,+rT +19- • iP (_,. ) 

Qm = Qm = 1 , Qm = -Qii^) = ^-^ — r 

nXkrl A{qk) 

(observe that from the L'Hopital-limit ^1^(1) = — ^^(1)), where 

n\) = 4(A)(A + (-l)'^^..,J(A-(-l)'^'.._,J, A(q) = a^-2a2cosq + l. (Ill) 

In the following we shall consider only the three cases which allow the normalization 
g^'^(O) = 1. Case (iii) can be treated too, but requires a special treatment, which here 
we shall not enter. According to which case the corresponding eigenvalue polynomial 
belongs, let us define the sets 'b^P\ V^p\ V^p^: 

k e V'^P^ if tP has a factor (A + r^Y, i.e. we have case (ii), 

k e V^P^ if tP has a factor (A - ruf, case (iii), and 

k e V'^p'> if tP has a factor (A^ — r^), i.e. we have case (iv). 

By D — \T>\ we denote the number of elements in D = U T>'^^\ similarly for V, etc. 
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7. Calculation of the matrix elements of in the homogeneous Ising model 

7.1. Explicit evaluation of the factors M{p') Ro{p') R{p') in (85) 

We now start to evaluate (85) with (83) and (86) for the homogenous Ising model where 
the parameters simplify drastically. Now 

C = h/a, rl = {a" - h'){a'- - l)/{a' - 1). (112) 

and u„ is represented by the Pauli cr^. 

We had agreed to consider initial states from the R-sector. Then for matrix elements 
of the final state must be NS. We specify the initial state by the momenta which are 
excited, i.e. by the 0"^ which are one, analogously the final state. Excluding for the time 
being case (iii), we take V^^^ to be empty. 

On the right of (83) we have to evaluate the factors A/'(p') Ro{p') R{p') ■ Let us 
start with i?o(p') = Qf\pi)Qf\pi) ■ 

For any choice of excitations, always one of the factors Qf\pi) or Qf\pi) is from 
case (i) of Subsection (6.4). Since we exclude for the moment case (iii), the other factor 
then must be from (ii) or (iv). So always Qf\Qi)Qf\Qi) — 1. For I e "D, case (n), we 
have Qf\l)Qf\l) = since either Qf\l) = or Qp^(l) = depending on the 
parity of /. So, in (83) the summation reduces to the summation over pi for Z e X> only, 
with fixed pi = for / G V. 

Ro{p') receives non-trivial contributions from Qk{^) of cases (i) and (iv). However, 
these can be written in a simple way if we use the explicit formulae for t*-^^(— r^) . For 
both values pi = 0, 1 the result is 

Q\%,) Ql"M = (-i)<"-"" tl)^ . Yl i-^rn + n ^^^^^ 

fi'T'i^i „ ri -\- T}^ 

where we get different results according to whether sq or or both (105) are excited: 



(-1) 
(-1) 



a2 - e'l 
62 e''^ - 1 



for (-1)^° = ±(-1)^- ; = (-l)^«<+l^l+'5i. (114) 



a2 - e'l 

Now, multiplying by Af{p') , it is easy to see that the products k & V in (113) cancel 
(recall that = for k eT) ) and we get finally 

1&T> m€V,m>l 

In the calculation of R{p') in (86) we have to insert our explicit expressions for 
t^'^^^Cn) ^iid t^^\(n) from (103). Here, as already mentioned after (105), the cases of 
even n and odd n give different formulas. E.g. the factor (A — (—1)°"'" s^r) is present only 
for R n even and NS n odd. So 

NS, nodd: t^'^-Q = (a'^ + !){-(+ (-irs^) J] ("C +rfc)^ H - rf) ,{116) 
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(for even n omit the bracket with s^), since in the NS-sector only odd k appear, and 
these fall into one of the classes (ii) and (iv), class (iii) being momentarily excluded. 
Analogously: 

R, nodd: t«(C) = (a2"-l)(C +(-ir«o) H (C + H -rf). (117) 

By slight manipulation we can move the pi-dependent terms such that they appear only 
in one place each in the numerator and get 

i?(-^^)(p') = R . {(-l)-(a^ + l)(-C +(-ir«.) Uievii-^rri +C) 

-(-ir(a2-i)(c + (-ir5o) a,p((-irn -0} ■ UkevH-^n +rk) im 

with 

R = (a/3)-("-'^/'a'^-^(a^" - l)/(a^ - 1), a = a' - 6^ /3 = 1 - a'fel (119) 

The first term in the curly bracket comes from the NS-sector final state, the second from 
the R initial state. The formula for Ri^^^^^) is similar. 

7.2. Summation, square of the matrix element 

Combining (115) with (118) the spin matrix element is given by a multiple sum over 
the components of p': 

($0kn|^l) 



( ^^0 I *0 ) 



E ( ^+ Uii-^y + + ^- n((-i) 

\ iGV lev 



pi 



n - C) X 



""11^^) n + Ci (-i)«n + (-i)p'«r^ ' ^^^"^ 

lev ' ^' meV,m>l V / * ' V / ^ 

with some p'-indcpcndcnt factors TZ'^. The superscript u is there to remind us that we 
have different expressions for n even and n odd, respectively. Now, in [27] it is shown 
that this sum can be performed, resulting in a factorized expression. As an example 
here we quote the summation formula for the multiple summation over pi with Z e D if 
the dimension of T> is odd and defined by the upper formula of (114): 

V- Uiev (-1)^' ((-i)^'n + ((-i)^'n + 

^fj^^ ni<m,i,mev ( (-1)'''^ + i-iy-r^) 

^ y (±(e^^'+'«'" - 1)) 

l,meV, l<m 
C = ct~^D-l){D-3)/4 ^-(I3-l)2/4_ 

The case of T> even is similar, see (52) of [27]. 
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In the following, we shall be interested in the product of the matrix elements of the 
spin operator between arbitrary periodic states, which does not depend on normalization 
of the left and right cigcnstatcs, i.e. we want to calculate 

(122) 





|u„| 




|u„| 






|<^'o)(<^'i 


l-^l) 



Taking the absolute squares, several factors in (121) can be rc-writtcn, e.g. 

1^2 gi^-, _ l|2 = |e'9' -a^p = A{qi) = - 2a^cosqi +1, (123) 

I e^^''+'^- - 1 P = ^^-^ A{q^) A{q,) '^'fj^^r\ (124) 
"P sm^ili - Qm) 

and all factors oc, f3 and A{q^) cancel. So we get for arbitrary n and ao — a^^: 





(7^ 

I n I 




< 




{ *o,o 


^0,0 







X 



,^ r ^-n ~ ^rn Sm + qm) 

Km, l,mGT> ^ ^ ' 

where 6 = 1. In a similar way we can find the product of matrix elements in the case 
of (To 7^ cr^. The final result is (125) with 5 = 0. Observe that the explicit appearance 
of excitations of type (ii), i.e. k & V has disappeared from our formula (recall that we 
still exclude k eV). 

7.3. Normalization of the periodic states, final result in terms of Aq, Xn, fk o-nd qi 

In order to compare (125) to the results obtained by A. Bugrij and O. Lisovyy [17, 18] 
we change the normalization and calculate the ratio (122). To do this we have to divide 
(125) by 

($o|<fo)(<fi|$i)/(^o,o|^o,o)'. (126) 

However, formula (76) cannot be used directly in our degenerate Ising case (102). As in 
the case (iv) we have first to go off the Ising point and consider ad — be = rj and apply 
rHopital's rule for — > 0. For n odd the result is 



(*o| 




l^l) 


(^o,o| 


1^0,0 


to 



k=l 



rifc-evenC'^Tr + rn,k) nfe-odd(''^0 + rn,k) 



X 



Ilk<i,k,i-odd{i^n,k + r„,z)(±r„,fc ± r„,,)) Ylk<i,k,i-even{K,k + r„,i)(±r„,fe ± r„,,)) 



and similar for n even, see 

Including also the hitherto excluded case (iii), our final formula for the matrix 
element is 

'r; + r„ s\n\{qi-q^) 





\<\ 


l^iX^il 


K\ 


l^o) 


($o|^o)(^i 


l^l) 



Kr. 
l,meT> 
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A„ Ukodd,ie.en(i-rk+ri){+rk-ri)) 



(127) 



where 



An = :. 7\^ TTo ^ ■ A f \ ^ TT-o ^ , odd 71, 



Hfegc^^) rifeeod) (^0 - 4) Ylkev^o) (A^+rfe) IlfeepW - n, 

nfcgp(°)(^+^fc)(A7r+rfc) 

" " ( Ao + A.) n,,^(i) ( Ao+r.) ( A,+r,) n,e^(^) ( - ( - ' 

Here we used a superimposed dot: ir^ as the short notation for if m e "D, for ±r^ 
if m e "D and for — if m e "D, respectively. For composite sets of momentum levels k 



we 



write V ^VUV, W = p(o)uP(°), ' = P«UP«. 



7.^. Fma/ result in terms of momenta 

Let {qi, q2, . . . , qx} and {pi, p2, . . ., pi,} be the sets of the momenta of the excitations 
presenting the states |$o) from the NS-sector and |$i) from the R-sector, respectively. 
After some lengthy but straightforward transformations of (127) we obtain 



(^0 














l^l) 



NS R 
K pNST-r2 AT L pRTT2 AT prif prL 

TT ^qfc -^Wlqfcl ^^q-qfe TT -^Pi lipT^IPil ^^P'Pi llfc=ill/=i^^qfc,Pi ^28) 



X 

fc=l 



where NS/2 (R/2) is the subset of quasi-momenta from NS (R) taking values in the 
segment < q < tt , NS/2 (R/2) containing q^; with odd k (even A;): 

_Sa±si sin^ si {si - si) 

"^"'^ " «a - sin ^ ' - {si - si) {si - si) ' 

NS NS R 

^ _ + g/3 _ riq^ (gp + gq) riq^ rip^gq + ^p)' 

^^Ce,l3 ) i7 R ■ NS R 

HpM^o+^p) nq!q,(sq + v)npV(^P+v) 

Fornodd: pNS fq ^ ^ _^ pR ^ fp ^ p _^ 

" (S^ - Sq)(So + Sq)' ' P (S^ + -Sp)(So - Sp)' 

TT — 

pR _ pNS _ 1 T _ lip' (-S^ + -Sp) ^ 

-^0 ~ -^TT ~ I ' ~ „ NS , , ' 

for n even: pNS ^ ^q pR ^ fp p ^ tt 

" (Stt + Sq)(SO + Sq)' ^ {Sn - Sp){So - Sp)' 



NS 



pR_ pNs_ 1 _ rip- + ■^,) 
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7. 5. Bugrij-Lisovyy formula for the matrix elements 

In [18] the following formula for the square of the matrix element of spin operator for 
the finite-size Ising model was conjectured: 



Ns( qi, q2, • • • , qi^ I < I pi, P2, • • • , pl ) 



Rl 



X 



n llq^qfcSmn ^ „ Ilp^p^siiiJi 2 h-T^y \ 

K -2 qfc-qfc/ L . 7 Pi-Pi' „:„1,2 7(qfc)+7(Pi) 

n ^^^ TT sm ^— -j-j- smh ^ 

-\nh^ 7(qfe)+7(qfcO 11 l{Pl)+l{Pi>) 11 sin^ ' ^ ^ 

fe<A;' smn 2 ^mn 2 i<fc<if 2 



In this formula the states are labelled by the momenta of the excitations. The factors 
in front of the right hand side of (129) are defined by 

where 7(q) is the energy of the excitation with quasi- momentum q: 

cosh 7(q) = _ , . , COS q , (130) 

and tx = tanhfC^, ty = tanhi^^^. 

Formula (129) can be easily derived from (128) if one takes into account the 
identification of parameters (14). In particular we have tx = ab, ty = {a — b)/{a + b) 
and the relation 

g7(q) ^ (131) 

a Sq — 

between the energy 7(q) of the excitation with quasi-momentum q and the 
corresponding zero Sq of the t(A)-eigenvalue polynomial (103). The following formulas 
give the correspondence between the different parts of (129) and (128): 

cc /f ".mh^ 7(q)+7(/?) f _f-l 

sinhi(7(0)+7(7r)) U-*x-V ' sin^ ^ tx-t'^ '^'^ " 

nNS ■ 7 qfc)+7 q , pNST-f— jvr 

q^q^ Smh 2 _ So + S^ ^q^ 1 lq^|q,| ^Vq,q, 



n sinh WpW sinh I^+IW jy^^^ 
R r^inVi 7(Pi)+7(p) I pR-rrf A 

p_^p^SmJl 2 _ So + Stt ^p, llp^lp^l ^V,Pi 



i.qfc 

nR r^inVi 7(Pi)+7(p) I pR-rrf ]\T 

P^Pl ^^^^ 2 ^ gp + S^ -^p, llp^lp,! ^^t 

For more details, see [27]. 
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8. Matrix elements for the diagonal-to-diagonal transfer-matrix and for the 
quantum Ising chain in a transverse field 

In this section we derive the matrix elements of the spin operator between eigenvectors 
of the diagonal-to-diagonal transfer-matrix for the Ising model on a square lattice (see 
Sect. 2.2). In this case the parameters are given by (15). As has been explained there, if 
we vary the parameters a and h in such a way to have fixed (a^ — 6^)/(l — a^6^) = 1/k', 
the eigenvectors (and therefore matrix elements) will not change. So we fix a = c = 
k'~^^^ and h — d — Q. Expanding the transfer-matrix (3) with such parameters we get: 



2A ^ ^ \ ^ 



2^ l^fc^fe+i + 9(^k): 
3 ^ k=i 

where Ti. is the Hamiltonian of the periodic quantum Ising chain in a transverse field. 
Prom (103) we get the spectrum of this Hamiltonian: 

£: = -ij]±£(q) (132) 
q 

where the energies of the quasi-particle excitations are 

£(q) = (l -2k'cosq + k'')^/2 ^ (|(k' - 1)^ + 4 k'sin' 1)'^' , qv^O, tt, 
e(0) = k'- 1, £{Tr) = k' + 1. 

In (132), the sign +/— in the front of e(q) corresponds to the absence/presence of 
the excitation with the momentum q. The NS-sector includes the states with an even 
number of excitations, the R-sector those with an odd number of excitations. The 
momentum q runs over the same set as in (103). Since we have a = c and b = d, the 
formula (128) with Sq = k'/£(q) for matrix elements for can be apphed. After some 
simplification we get the analogue of (129), now for the quantum Ising chain: 

|Ns(qi,q2,...,qi.K|pi,P2,...,PL)R| =k' ^ Ji ji X 

where 




.=1 1=1 V 2 sin 2 



nf np(^(q)+^(p))^ 



n^^' (^(q) + ^(q'))^ nSp' (^(p) + ^(p'))^ 



and 



nj'(£(q)+£(q')) 



e 



n^(6(q)+6(p)) ■ 

Formally, all these formulas are correct for the paramagnetic phase where k' > 1, and for 
the ferromagnetic phase where < k' < 1. But for the case < k' < 1 it is natural to 
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redefine the energy of zero-momentum excitation as £(0) = 1 — k' to be positive. From 
(132), this change of the sign of £(0) in the ferromagnetic phase leads to a formal change 
between absence-presence of zero-momentum excitation in the labelhng of eigenstatcs. 
Therefore the number of the excitations in each sector (NS and R) becomes even. Direct 
calculation shows that the change of the sign of £(0) in (133) can be absorbed to obtain 
formally the same formula (133), but with new £(0) and even L (the number of the 
excitations in the R-sector) and new ^ = (1 — k'^)^/^. 

Formulas (129) and (133) allow to re-obtain well-known formulas for the Ising 
model, e.g. the spontaneous magnetization [1, 2]. Indeed, for the quantum Ising 
chain in the ferromagnetic phase (0 < k' < 1) and in the thermodynamic limit 
n ^ oo (when the energies of |vac)Ns and |vac)R coincide, giving the degeneration 
of the ground state), we have — 1 and therefore the spontaneous magnetization 
Ns(vacK|vac)R = = (1 - k'^)V8. 

9. Conclusions 

We have shown that finite-size state vectors of the Ising (and generalized Ising) model 
can be obtained using the method of Separation of Variables and solving explicitly 
Baxter equations. The Ising model is treated as a special N = 2 case of the Zjy-Baxtcr- 
Bazhanov-Stroganov r*^^^-model. Finite-size spin matrix elements between arbitrary 
states arc calculated by sandwiching the operators between the explicit form of the 
state vectors. For the standard Ising case this gives a proof of the fully factorized 
formula for the form factors (129) conjectured previously by Bugrij and Lisovyy. 

We also extend this result to obtain a factorized formula for the matrix elements 
of the finite-size Ising quantum chain in a transverse field. We show how specific local 
operators can be expressed in terms of global elements of the monodromy matrix. The 
truncated functional relation guaranteeing non-trivial solutions of the Baxter equation 
is compared to Baxter's [43] Ising model functional relation. 
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